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Abstract— Decentralized power generation is an option to serve 
communities that survive in places without electricity or isolated from 
urban centers or in areas with intermittent electricity coming from gasoline 
or diesel generator sets. This article presents an alternative energy 
generation with low environmental and financial impact, a floating 
hydroelectric generation system using an axial flow electric generator 
combined with a cross-flow hydrokinetic converter that works from water 
flows in rivers’ sections floating structure allows for performance even at 





I. INTRODUCTION 


Electricity is needed to develop isolated regions of the 
country, as access to this service is vital in economic and 
social development, impacting people's communication, 
productive use, and a direct relationship to improving the 
quality of life [4 ]. 


The high costs of transmission lines can be caused by 
the large distances between homes to the network due to 
the low population density in the Amazon region [4]. 


Alternative decentralized service solutions have been 
encouraged by the federal government through social 
programs. These programs aim to bring electricity to rural 
regions, isolated through micro and mini hydroelectric 
plants, thermoelectric plants powered by biofuels or 
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natural gas, solar photovoltaic plants, wind turbines, 
hybrid systems [3], [4]. 


The proposed Hydrokinetic Power Generation System 
coupled to the Axial-Flow Generator may be a viable 
option for communities close to rivers with specific 
hydrological characteristics, including a frequent and 
constant water flow. This generation system is robust, easy 
to maintain, and with low financial investment and no 
environmental impact, as it does not use large, flooded 
areas. 
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II. MODELING OF THE HYDROKINETIC 
GENERATION SYSTEM COUPLED TO THE 
AXIAL FLOW GENERATOR 


Michell-Banki-type cross-flow turbines use a drum- 
type rotor equipped with blades. The water flow is directed 
by a hydraulic profile and meets the rotor blades, driving 
it. This type of turbine is used in hydraulic projects with 
smaller waterfalls and a greater volume of excellent water, 
which usually occurs in flatter regions. [1], [3]. 


The various types of couplings used between the 
converter and the generator result in problems ranging 
from balancing the axes, exposure of the parts to weather, 
more significant space used, and even maintenance [3]. 


The construction of unique equipment, which combines 
a cross-flow hydro converter with a permanent magnet 
generator or squirrel cage induction motor in an axial flow 
machine, allows for space and maintenance gains and 
greater power [3]. The proposed system works submerged 
or in line with pipes where there are flows of any type [3], 
[5]. 

In [3] presents the mathematical modeling of the axial 
flow generator and the dimensioning of the permanent 
magnets. 


The prototype will be installed on the Madeira River. 
One of the largest rivers in the world in terms of water 
volume. It has a total drainage area of 1,420,000 km2, or 
about 1/5 of the entire Amazon basin. 


The surveys were carried out with the National Water 
Agency (known in Brazil by its Portuguese acronym - 
ANA), in its database, and with the Mineral Resources 
Research Company (known in Brazil by its Portuguese 
acronym - CPRM) in Rondônia, where there is historical 
data on flow, flow velocity, and depth, were made 
available for consultation, enabling modeling and 
simulation for the elaboration of the construction project of 
the floating hydroelectric system in realist scale. 


Fig.l shows the installation location of the 
Hydrokinetic Power Generation System coupled to the 
floating Axial Flow Generator. 





Fig.1: The Installation location of the proposed system. 
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The Madeira River is known for transporting a large 
amount of suspended material and, in the rainy season, its 
flow goes from 5,000 m°/s to 50,000 m?*/s. Understanding 
the dynamics of the river, given this seasonality, is 
essential for determining the parameters of the floating 
hydroelectric system. 


A large amount of sediments in the river increases the 
fluid density and, in a direct ratio, the available hydraulic 
power. However, the corrosion problems are more evident. 


An important input parameter for the design of a 
hydrokinetic converter is the river flow velocity since the 
available hydraulic power is directly proportional to the 
section area and the cube of the flow velocity. This data 
over time becomes even more relevant in a river with 
significant variation in its regime, as is the case of the 
Madeira River, as shown in Fig. 2. 


PORTO VELHO 
Fluviermetrica - 15400000 
Vazao- 18 





Fig.2: In red is a historical maximum, in orange historical 
minimum, and in blue the flow value of the last 12 months. 


An ACDP - Acoustic Doppler Current Profiler is used 
to raise the current profile, a hydro-acoustic current 
profiler that measures the velocity of particles in the water 
column using a physical principle of propagation of sound 
waves (Doppler effect). This type of equipment is shown 
in Fig. 3. 





Slow+— Fast 








Fig.3: Functioning Acoustic Doppler Current Profiler [6]. 


The collected flow velocity historical data are 
simulated in the WinRiver II software, where it is possible 
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to observe the current velocity profile and the bed trace as 
shown in Figures. 4 and 5. 
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Fig.4: Typical measurement in the flood period (February 
2020). 
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Fig.5: Typical measurement in the dry period (September 
2016). 


As shown in figures 4 and 5, the average velocity at 
200 meters from the right bank of the Madeira River varies 
from 0.7 m / s in the dry period to 2.5 m / s in flood. With 
this data, it is possible to develop simulations of the river 
flow, including the floating system where the turbine will 
be located. Fig. 6 shows the flow simulation of a natural 
flow channel with the floating turbine housing system. 

Pwnage TTO- 00T- 





Fig.6: One-channel flow simulations. 


A simulation of a diffuser to determine the energy 
available for a generation was presented in [3]. The results 
presented demonstrate an extensive recirculation below the 
diffuser. This recirculation is a hydraulic loss. A new 
diffuser geometry is proposed to solve this problem, 
according to the simulations presented in Figures 7, 8, 9, 
and 10. 
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Fig.7: Front view of the diffuser. 





Fig.8: Side view of the diffuser. 
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Fig.9: Rear view of the diffuser with exit window (rotor 
position). 





Fig.10: Isometric view of the diffuser. 


The computational model used for simulation is shown 
in Fig. 11. The inflow (upstream), Fig. 12, was based on 
historical measurements of the speed of the Madeira River 
in Porto Velho — Rondônia. The geometry was worked on 
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the SpaceClaim software, Ansys Inc. 2017, and the 
Meshing software generated the mesh. 


In the dry period, the average speed is around 0.7 m/s, 
and in the flood period, it is 2.5 m/s. Based on this history, 
the input speed was simulated from 0.5 to 3.0 m/s with an 
increment of 0.5 m/s. 





000 150) 1,008 (mn) 
a 





Fig.11: Computational model used for simulations. 
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Fig.12: Computational model boundary conditions. 


For numerical simulation, used the CFX software 
Ansys Inc. 2017. The parallel processing resource was 
used. 


In Fig.12, the input boundary conditions (upstream) are 
presented along with the average speed obtained near the 
diffuser inlet and the average diffuser output speed. 


In Table 1, the input boundary conditions (upstream) 
are presented along with the average velocity obtained 
near the diffuser inlet and the average diffuser exit 
velocity. 


Table I - Speeds Obtained in the Computational Model 


Input Average input speed Average output speed 


3,0 2,880090 2,998580 
2D 2,400370 2,497730 
2,0 1,920250 1,997130 
1,5 1,440110 1,496750 
1,0 0,959991 0,996668 
0,5 0,479909 0,497159 
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The input velocity is the boundary condition imposed 
on the computational model, which, in turn, represents the 
historical average velocity measured in the Madeira River. 
Due to the flow pressure drop and viscous and kinetic 
effects, the speed decreases along the way. An average 
speed near the diffuser inlet (blue plane) was calculated, 
Fig. 13, which shows the details of the regions to obtain 
the average velocities. 


The hydraulic potential across the diffuser is obtained 
by an average speed calculated at the outlet (yellow plane). 
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Fig.13: Details of regions to obtain average speeds. 
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Fig.15: Velocity vector field detail for 0.5 m/s flow. 


Figures 14, 15, and 16 show a significant reduction in 
recirculation below the diffuser. 
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Fig.16: Velocity field for 0.5 m/s flow. 


As can be seen, there is a high-velocity region below 
the diffuser. By the principle of mass conservation, the 
tendency is for the flow to pass through regions with less 
obstruction. Nevertheless, we observed a loss of speed 
along the path (Fig.16); however, we observed that the 
diffuser played its role, increasing the rate at its exit. The 
figures below show different velocity profiles in the 
equipment as a function of different initial flow velocities. of 
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Fig.18: Velocity vectors for input flow of 1.0 m/s. 
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Fig.22: Velocity field for 1.5 m/s flow. 
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Fig.24: Velocity vector field detail for 2.0 m/s input flow. Fig.28: Velocity field for 2.5 m/s flow. 
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Fig.26: Velocity vectors for input flow 2.5 m/s. Fig.30: Velocity vector field detail for 3.0 m/s input flow. 
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Fig.31: Velocity field for 3.0 m/s flow. 


The model presented will be printed (3D printer) in 
reduced scale for testing in the Control and Automation 
laboratory channel to validate numerical results. As only 
the flow of the Madeira River is known, through the 
simulation, an average rotation that would be produced by 
the flow passing through the diffuser was estimated. We 
will use this rotation to simulate both the electrical and 
mechanical parts of the rotor to obtain the shaft torque 
produced in the turbine. 


It determined an average rotation of the average speed 
at the diffuser output. Considering that the rotor blade is 1 
m long and the linear velocity is proportional to the 
angular velocity: 


V = wR (1.0) 


Table 2 presents the values of the theoretical average 
rotation obtained from the average output speed of the 
rotor. 


Table 2: Theoretical average rotation obtained at the 
diffuser output 








Input Average Rotation (r.p.m.) 


3,0 57,2686595 1 
2,5 47,10312912 
2,0 38,14237338 
1,5 28,58581933 
1,0 19,03495666 
0,5 9,495037482 








This average rotation will be used in the modeling of 
the electrical part and the simulation of the rotating part of 
the turbine (rotor). 


The generator's electrical machine can either be built 
with permanent magnets or a squirrel-cage induction 
motor, and the voltage generated varies according to the 
installation's water flow. A proprietary conversion system, 
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built into the system body, synthesizes frequency and 
voltage to desired values. This generation system will be 
made with permanent magnets. 


Considering the movement of the permanent magnet 
rotor, the distance traveled by a pair of magnets induces a 
voltage cycle at the terminals of a stator coil taken as a 
reference. Thus, adopting thirty magnets on the concentric 
segment of magnetizable material, adjacent to the face 
with fixed magnets, spaced by a few millimeters, there will 
be another concentric segment for the mechanical support 
of the nine stator units, generating different phases of 
electrical voltage. 


Due to the initial perspective of adopting three 
concentric stators in the full-scale design, there would be 
different tangential speeds for the same rotor speed. A fact 
that inevitably conditions the design of the coils to three 
distinct types. In the smaller-scale design, each of these 
three stator coil designs is on a single circle. There should 
be different assessment criteria, adopting three speeds for 
each test. 


The rotation speeds in the reduced-scale prototype will 
be such that they must induce voltage at the frequencies of 
35.2 Hz, 29.4 Hz, and 25.2 Hz, respectively, those 
observed in the full-scale model design, from the most 
external stator to the more internal. 


Di C m D >C S BASE ole 








ome. Z gt 
R ‘ F/F RN “tna 








Fig.32: Small-scale prototype stator and rotor. 


Considering that the passage of a pair of magnets from 
different magnetic poles induces an electrical voltage cycle 
at the unit stator terminals taken as a reference, 
considering the rotor with 30 magnets, for those steady- 
state frequencies, follows the proportion below: 


(35,2 Hz’ 
(293 z) (X) magnets 
25,2 Hz. 20 
1 Hz 2 magnets (An 
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Therefore: 


704 
X= (503) magnets /s 


50.4 (3.0) 


Thus, the steady-state will be reached from the 
outermost to the innermost stator in the full-scale model. 
There are displacement rates of 70.4, 58.8, and 50.4 
magnets per second. 


Considering the model at a minimum scale with one 
circular path that containing the thirty magnets, the 
necessary velocities for each of the three different stator 
units to reach the steady-state will be obtained with the 
proportion below: 


‘70,4 
(503) magnets js 
50,4 


1 rev iP 


2,35 "EP /, ===> = 141,0 RPM 





30 magnets /s 1,96 T /s = — = 117,6 RPM 
168 TEY /s = => = 100,8 RPM 


(4.0) 


The rotations calculated for the model in minimum 
scale are intended to validate the theoretical and practical 
arrangement, leading to the refinement of the electrical 
design of the machine in full scale. 


For the best understanding of the study, each of the 
three stator unit designs will be fixed side by side in a 
single path of excitation by the rotor. They will be 
identified as R-phase, S-phase, and T-phase. The values 
for phases R, S, and T represent the stator rotation on a 
minimum scale with the same speed of passage of the 
magnets provided in the full-scale model. 


Phase R 


e Rotation in steady-state 141.0 r.p.m.; 

e Number of turns per stator unit 85 (15 AWG); 

e Maximum electrical current in conductors 
(450 A/cm?) 7,427 A; 

e Theoretical value (per stator unit) at 16.0 V 
steady state rotation. 


Phase S 


e Rotation in permanent regime 117,6 r.p.m.; 

e Number of turns per stator unit 109 (15/16 
AWG); 

e Maximum electrical current in conductors 
(450 A/cm’) 5,890 A ; 

e Theoretical value (per stator unit) at 17.2 V 
steady state rotation. 

Phase T 


e Rotation in steady-state 100.8 r.p.m.; 
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e Number of turns per stator unit 136 
(17 AWG); 

e Maximum electrical current of conductors 
(450 A/cm’) 4,671 A; 

e Theoretical value (per stator unit) at 18.5 V 


steady-state rotation. 


IHI. PROTOTYPE OF THE HYDROKINETIC 
GENERATION SYSTEM COUPLED TO THE 
AXIAL FLOW GENERATOR 


The model was printed on a reduced-scale 3D printer 
for testing in the channel belonging to the Control and 
Automation laboratory to validate the results presented in 
the simulation. 





Fig.33: Hydrokinetic Power Generation System coupled to 
the Axial-Flow Generator. 





Fig.34: Turbine and diffuser overview. 
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10:09:25 
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Fig.35: The three-phase voltage generated in loads of I Q 
of 30 W. 


IV. CONCLUSION 


The new diffuser geometry proposed in this article 
increased its output speed as predicted in the simulations 
and the test channel. 


The built prototype can generate energy for isolated 
communities and can be used in various economic 
activities with rivers, canals, or small watercourses, with 
low flow speeds. 


The proposed system is relevant for its low 
manufacturing and installation cost. Furthermore, it 
contributes to the preservation of the environment, as there 
is no need for dams. 
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